Anisotropic optical response of the mixed-valent Mott-Hubbard insulator NaCu 2 C>2 



O 

(N 



•4— > 

a 

C 

o 
o 



> 
o 

(N 



On 
O 



Y. Matiks, 1 ^ A. N. Yaresko, 1 K. Myung-Whun, 1 - 2 A. Maljuk, 3 ' 1 P. Horsch, 1 B. Kcimer, 1 and A. V. Boris 1 '!! 

1 Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany 
2 Department of Physics and IPIT, Chonbuk National University, Jeonju 561-756, Korea 
'Leibniz Institute for Solid State and Materials Research IFW Dresden, Helmholtzstr. 20, D-01171 Dresden, Germany 

(Dated: September 20, 2011) 

We report the results of a comprehensive spectroscopic ellipsometry study of NaCu2 02, a com- 
pound composed of chains of edge-sharing Cu 2+ 04 plaquettes and planes of Cu 1+ ions in a 0-Cu 1+ -0 
dumbbell configuration, in the spectral range 0.75 — 6.5 eV at temperatures 7 — 300 K. The spectra 
of the dielectric function for light polarized parallel to the Cu 1+ planes reveal a strong in-plane 
anisotropy of the interband excitations. Strong and sharp absorption bands peaked at 3.45 eV (3.7 
eV) dominate the spectra for polarization along (perpendicular) to the Cu 2+ 02 chains. They are 
superimposed on flat and featureless plateaux above the absorption edges at 2.25 eV (2.5 eV). Based 
on density-functional calculations, the anomalous absorption peaks can be assigned to transitions 
between bands formed by Cu 1+ 3d xz (d yz ) and Cu 2+ 3d xy orbitals, strongly hybridized with O p- 
states. The major contribution to the background response comes from transitions between Cu + 
Sd z 2 and 4p x (p y ) bands. This assignment accounts for the measured in- plane anisotropy. The dielec- 
tric response along the Cu 2+ C>2 chains develops a weak two-peak structure centered at 2.1 and 2.65 
eV upon cooling below ~ 100 K, along with the appearance of spin correlations along the Cu 2+ 02 
chains. These features bear a striking resemblance to those observed in the single- valent Cu 2+ 02 
chain compound LiCuVCU, which were identified as an exciton doublet associated with transitions 
to the upper Hubbard band that emerges as a consequence of the long-range Coulomb interaction 
between electrons on neighboring Cu 2+ sites along the chains. An analysis of the spectral weights of 
these features yields the parameters characterizing the on-site and long-range Coulomb interactions. 
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I. INTRODUCTION 



Copper oxides with quasi-one-dimensional electronic 
structure have drawn much attention because of their 
unusual magnetic properties and the variety of ground 
states originating from strong electronic correlations in 
different lattice architectures. Compounds composed of 
chains of edge-sharing Cu 2+ C>4 plaquettes have the pecu- 
liar property that the magnitude of the nearest-neighbor 
hopping matrix clement along the chains is anomalously 
small due to the orthogonality of the 2p a orbitals on 
oxygen ions adjacent to the Cu ion. The interplay be- 
tween short- and long-range interactions generates spiral 
magnetism in Mott-insulating compounds*— and charge 
density waves in doped compounds4£ By virtue of their 
exceptionally narrow electronic bandwidths, these com- 
pounds also provide a highly favorable platform for the 
investigation of exciton formation and the interplay be- 
tween spin and charge correlations in the cuprates£~— 

In optical experiments, Zhang-Rice singlet excitations 
are generated in the charge transfer process d^d^ —> 
(d 9 Lh)idj° between Cu 2+ 04 plaquettes at sites i and 
jiii In the final state, the spins form singlets leaving an 
oxygen- ligand hole on one of the two Cu 2+ 04 plaquettes. 
The band formed by the d 10 states on the other plaque- 
tte can be regarded as the upper Hubbard band, in for- 
mal analogy to the single-band Hubbard model. Using 
spectroscopic ellipsometry, Matiks et al. have recently 
demonstrated that Mott-Hubbard excitons are the low- 
est accessible states for holes in the single-valent chain 



cupratc LiCuVC^ii 2 - For photon polarization along the 
Cu 2+ 02 chains, a weak but well-resolved two-peak struc- 
ture whose spectral weight is strongly enhanced upon 
cooling near the magnetic ordering temperature has been 
identified as an exciton doublet. These results have not 
only persuasively demonstrated the formation of Mott- 
Hubbard excitons, but also quantified various impor- 
tant characteristic energy scales, such as the local Hub- 
bard and long-range Coulomb interactions, the nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) hop- 
ping parameters, and the resulting superexchange ener- 
gies. 

In order to explore the generality of the excitonic states 
observed in LiCuV04, we have carried out a comprehen- 
sive ellipsometric study of the complex dielectric func- 
tion of single crystals of NaCu202, a mixed-valent Mott 
insulator composed of chains of edge-sharing Cu 2+ 04 
plaquettes and 0-Cu 1+ -0 dumbbells*! 3 . Like LiCuV04, 
NaCri202 exhibits helical magnetic order at low temper- 
atures due to competing ferromagnetic NN and antifer- 
romagnetic NNN superexchange interactions within the 
Cu 2+ 02 edge-sharing chainsJ^r— Unlike in LiCuV04, 
the chains in NaCu202 are separated by nonmagnetic 
Cu 1+ ions in 0-Cu 1+ -0 dumbbell complexes along the 
c axis. The presence of this additional structural unit 
in other copper oxides including YBa2Cu306, the parent 
compound of a well known family of superconductors, 
provides further motivation for our study. 

In a recent ellipsometry study of LiCu2C>2, a compound 
that is isostructural and isoelectronic to NaCu202, Pis- 
arev et al. reported a strong and narrow absorption 
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peak at about 3.3 eV and attributed this excitation to 
the 0-Cu 1+ -0 dumbbellsiii This feature obscures the 
weak exciton bands in the optical conductivity along 
the Cu 2+ C>2 chains that are expected following the anal- 
ogy with LiCuV04. Detailed information about the 
in-plane dielectric anisotropy is required in order to 
clearly separate the contribution of the Cu 2+ (d 9 ) and 
Cu 1+ (d 10 ) states to the dielectric response and to avoid 
ambiguity in the data interpretation. However, Li-Cu 
chemical inter-substitution and crystallographic twinning 
in LiCu202 crystals obliterate the dielectric anisotropy 
in the ab plane. Experiments on single crystals of 
NaCu202 offer the chance to elucidate the intrinsic di- 
electric anisotropy of this class of compounds. Indeed, 
recent x-ray spectroscopy and neutron diffraction investi- 
gations confirmed the superior quality of NaCu2C>2 single 
crystals which, unlike LiCu2 02, are not prone to twin- 
ning and disorder J^"— The absence of crystal defects can 
explain the apparent absence of multiferroic behavior in 
NaCu 2 2 ^^ 

In this paper, we report a comprehensive ellipsometric 
study of the dielectric function anisotropy of NaCu2C>2 
in the spectral range 0.75 - 6.5 eV and its interpretation 
based on band-structure calculations. The paper is orga- 
nized as follows. Section II describes experimental details 
and results. In Section III A the discussion is focused on 
the behavior of the exciton bands identified in the optical 
conductivity along the Cu 2+ C>2 chains. The anomalous 
optical absorption is discussed in Section III B, and in 
Section III C local spin density approximation (LSDA) 
calculations are reported in order to explain the observed 
anomalies and the anisotropy of the optical response. Fi- 
nally, our conclusions are summarized in Section IV. 



II. EXPERIMENTAL DETAILS AND RESULTS 

Single crystals of NaCu2 02 were grown by the self- 
flux technique, as described in Ref. [lj. X-ray diffrac- 
tion and inductively coupled plasma atomic spectroscopy 
measurements showed no impurity phases and a chem- 
ical composition consistent with ideal stoichiometryj^ 
The crystal structure of NaCu202 belongs to the Pnma 
space group with an orthorhombic crystal structure and 
room-temperature lattice parameters a = 6.2087A, b = 
2.9343A, c = 13.0648AJ£ The unit cell is composed of 
two pairs of the edge-sharing Cu 2+ C>2 chains, running 
along the b axis and shifted relative to each other by 
6/2, see Fig. 1. There are four Cu 2+ ions per unit cell, 
and the density of Cu 2+ ions is N Cu 2+ = 1.68 x 10 22 
cm~ 3 . Within the single chain, the Cu 2+ — Cu 2+ distance 
is 2.934 A, and the Cu 2+ -0-Cu 2+ bond angle is 92.9°. 

Magnetic susceptibility data follow the Curie- Weiss 
law above 150 K, with a Curie temperature of —62 KJ^ 
A broad maximum typical for low-dimensional magnets 
is observed at T ~ 52 K, associated with short-range 
antifcrromagnetic correlations within the chains. Low- 
temperature susceptibility and specific heat data indicate 




Figure 1: Crystal structure of NaCu2C>2. 

a magnetic phase transition at = 11-5 Kf^ Because 
of the 92.9° Cu— O— Cu bond angle within the chains, 
the NN exchange integral is small (Ji = —16.4 K) and 
dominated by that between the next-nearest neighbors 
(J2 = 90 K). Based on a refinement of single-crystal neu- 
tron diffraction data, an elliptical helix structure with 
alternating polarization planes was proposed as the mag- 
netic ground state.— 

For our optical measurements a single crystal of 
NaCu202 with dimensions 5x5x0.2 mm and the freshly 
cleaved ab surface was used. The cleaving and mounting 
procedures were performed in argon atmosphere to pre- 
vent sample oxidation. The transfer of the cleaved and 
mounted sample into the cryostat was carried out in an 
argon-flooded glove bag. The sample was mounted on the 
cold finger of a helium-flow cryostat with a base pressure 
of 2 x 10 -9 Torr at room temperature. The ellipsometric 
measurements were performed with a rotating-analyzer 
type Woollam VASE variable angle cllipsomcter. 

Room temperature spectra of the ellipsometric angles 
^(oj) and A(oj) are shown in Fig. 2. The angles ^ and A 
are defined through the complex Frcsncl reflection coeffi- 
cients for light polarized parallel (r p ) and perpendicular 
(r s ) to the plane of incidence, 

tan #e iA = r p /r s . (1) 

The data were measured at angles of incidence, tpi, of 
65° and 72.5° for sample orientations with the a or & 
axis in the plane of incidence, which denoted as _L and 
||, respectively, corresponding to the orientation of the 
electric field with respect to the Cu 2+ 02 chains. The 
observed dependence of the ellipsometric spectra on the 
sample orientation points to an intrinsic anisotropy of the 
optical response of NaCu202- 

Ellipsometry yields the anisotropic frequency- 
dependent complex dielectric tensor, 

e(w) = si(uj) + ie2{uj) = l + i 4k<t{<jj)/uj, (2) 

without the need for reference measurements or Kramers- 
Kronig transformations. A numerical regression 
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Figure 2: Ellipsometric angles ty(u)) and A(ui) measured at 
room temperature for ab sample surface orientations with the 
(a) a or (b) b axis in the plane of incidence. 
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Figure 3: Room temperature spectra of the real ei and imagi- 
nary £2 parts of the dielectric function for polarizations along 
and perpendicular to the chains of NaCu202. For polariza- 
tion along the chains, the spectra are shifted by 2 for ei and 
4 for ei for clarity. Dotted lines depict the pseudo-dielectric 
functions. 



procedure^ was applied to derive the principal compo- 
nents of the dielectric tensor from the ellipsometric data 
at different Euler angles. Since measurements of the op- 
tical response along the c axis could not be obtained 
due to the thinness of the sample, spectra measured at 
two angles of incidence for each polarization were used in 
the fitting procedure. The result of fits to the data are 
shown by thin lines in Fig. 2. The real and imaginary 
parts of the dielectric function for both polarizations re- 
sulting from the fit are drawn in Fig. 3. In this figure 
we also show the pseudo-dielectric function which is ob- 
tained by direct inversion of Eq. (1) for >F(u;) and A(ui) 
measured at ipi = 65°, assuming isotropic media with 
£xx = £yy = £zz- The effect of anisotropy mainly ap- 
pears in the pronounced peaks in the spectral range 2.5 
- 4.0 cV. From this point on we will discuss only the 
corrected dielectric function. 

The optical spectra for both polarizations can be bro- 
ken up into a, /? and 7 zones, as shown in Fig. 3. The 
a zone comprises the spectral range of the strong and 
narrow absorption peaks located at 3.45 eV (along the 
chains) and 3.7 cV (perpendicular to the chains) that 
are dominating the spectra. The anomalous strength of 
these excitations leads to negative values of e\. The low- 
energy shoulders of these peaks near the insulating gap 
are marked as the ft zone. The 7 zone includes the broad 
and weak features above 4 eV that give rise to a rather 
flat and featureless optical response. 

The temperature dependencies of the real parts of the 
optical conductivity ui{uj) and the real parts of the di- 



electric function £i(w) for both polarizations are plotted 
in Fig. 4. The strongest temperature effect on the optical 
spectra observed is a growth of the sharp 3.45 (3.7) cV 
peak for polarization along (perpendicular to) the chains 
with decreasing temperature. Against the background of 
this growth, some changes are apparent on the low-energy 
shoulders of the strong peaks. The low-temperature fea- 
ture centered near 2.65 eV for polarization along the 
chains (see inset of Fig. 4(b)) is the most intriguing of 
these features. 

To study the temperature-driven changes in detail, 
temperature-difference spectra A<7i(w,T) and Ae\(ui,T) 
with respect to 150 K are plotted in Fig. 5. Perpen- 
dicular to the chains, ostensible changes in Ae^(u>) and 
Aa^(ui) are seen in the a and (3 zones between 2.5 and 4.5 
cV. Alternating areas of positive and negative regions in 
the Aa^-(u>) spectra within this spectral range are mutu- 
ally compensated. This observation indicates that the 
individual optical-band intensities within this spectral 
range are conserved, and that the spectral weight is not 
redistributed among these bands. Rather, the changes 
are induced by concurrent narrowing of the pronounced 
band at 3.7 eV and the optical band located on its low- 
energy shoulder at 3.37 cV upon cooling. Above 4.5 eV 
there are no other changes, except minor variations near 
5.5 eV that can be attributed to a slight narrowing of the 
5.4 eV band. 

Along the chains, the appearance of the feature marked 
with an arrow in Fig. 5(b) makes the temperature- 
difference spectra Ae\(u>) and Acr|(a;) qualitatively dif- 
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Figure 4: Real part ei of the dielectric function and real part 
el of the optical conductivity for polarizations along and per- 
pendicular to the chains of NaCu2C>2 measured at different 
temperatures. Inset: Magnified view of o\ and £i for polar- 
ization along the chains. 



ferent from those for polarization perpendicular to the 
chains. The peak in Aerf(u;) at 2.65 eV and anti- 
resonance with zero-crossing at the same energy in 
Aef (ui), which develop below 150 K, indicate an enhance- 
ment of the optical band upon cooling. The changes of 

the erf amplitude with temperature at the peak position 

for this band and ef right below and above, as marked 
with vertical dashed lines, were studied by temperature 
dynamic scans shown in Fig. 5(c). A Kramcrs-Kronig 

consistent enhancement of erf at 2.62 eV below 150 K, 

together with an upturn and a downturn in e\ at 2.45 
and 2.74 eV, respectively, confirm the pronounced inten- 
sity increase of the 2.65 eV band with decreasing tem- 
perature. 

The transfer of spectral weight can be quantified by 
integrating the optical conductivity in terms of the effec- 
tive charge density 
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Figure 5: Temperature-difference spectra Aeri(w,T) = 
ai(td,T)-(Ti(u), 150 K) and Asi(u,T) = ei(u,T) -si(td, 150 
K) of NaCu202 for polarizations (a) perpendicular to and (b) 
along the chains, (c) Temperature dependence of erf measured 
at 2.62 eV and ef measured at 2.45 and 2.74 eV for polar- 
ization along the chains, as marked by vertical dashed lines 
in (b). (d) Spectral weight changes AN e ff(id) = N e ff(uj,7 
K) — N e ff(id, 100 K) for polarization along the chains. 



where m is the free electron mass. Figure 5(d) shows 
the changes in the spectral weight below 100 K along the 
chains. The rise of the spectral weight at low energies 
is due to the band at 2.65 eV. Above 3 eV this growth 
is compensated by a loss in the spectral weight of the 
higher-energy bands, following the optical sum rule. The 
feature marked with a rectangle in Fig. 5(d) originates 
from the narrowing of the strong optical band at 3.45 
eV with decreasing temperature and indicates that the 
spectral weight is retained within this band. 

To separate the contributions from different optical 
bands to the dielectric function spectra for both polar- 
izations, a classical dispersion analysis was performed. 
A minimum set of Lorcntzian oscillators, with one high- 
energy oscillator beyond the investigated spectral range, 
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Figure 6: Real ei and imaginary £2 parts of the dielec- 
tric function of NaCu2 02 measured at 7K (black thick line) 
and dispersion analysis fit to the data (light thin line) for 
photon polarizations (a) perpendicular to and (b) along the 
Cu 2+ 02 chains. Thin color lines represent the principal 
bands contributing to the optical response (see Table I), 
as derived from the dispersion analysis, (c) Temperature- 
difference spectra Aei(u;) = £i(w, 7 K) — £i(oj,150 K) and 
Ae2(w) = £2(w,7 K) — £2(<^,150 K) for polarizations along 
the chains of NaCu202, and results of 'fit 1' and 'fit 2', as 
described in text, fnsets: (a) Temperature dependence of the 
bandwidths of the Lorentz oscillators at 3.45 eV (black) and 
3.7 eV (red) for polarizations along and perpendicular to the 
chains, respectively, (b) Temperature dependence of the os- 
cillator strength of the 2.65 eV band for polarization along 
the chains. 



was introduced to represent a dielectric function in the 
form 

e fro= e c,+x; , j j _ , (4) 



where Wj , Tj , and Sj are the peak energy, width, and os- 
cillator strength of the jth oscillator, and £oo is the high 
energy core contribution to the dielectric function. The 
parameters of the individual complex Lorentzian oscilla- 
tors were derived with high accuracy and reliability by 
simultaneously fitting to £i(w) and e^ifjj)- They are listed 
in Table I. 

Figures 6(a) and 6(b) represent the principal bands 
contributing to the dielectric response at 7 K for polariza- 
tions along and perpendicular to the chains. According 
to the dispersion analysis, the a zone of the spectra is 
composed of two narrow bands at 3.63 (3.40) and 3.71 
(3.46) eV for polarization perpendicular to (along) the 
chains, giving in sum the extremely strong asymmetric 
absorption band at 3.7 (3.45) eV. The 7 zone of the spec- 
tra is composed of a series of excitation bands. Along 
the chains, four nearly equally spaced optical bands (at 
4.01, 4.41, 4.75 and 5.16 eV) can be recognized. Per- 
pendicular to the chains, the dispersion analysis gives a 
broad band at 4.85 eV. However, this band can be split 
into two subbands that are surrounded by excitations at 
4.30 and 5.40 eV, forming a series of excitations similar 
to that for polarization along the chains. The 7 regions 
in both polarizations end in a broad weak band near 6 
eV and one optical band beyond the investigated spectral 
range. Only one band located at 3.37 eV is lying in the f3 
zone for polarization perpendicular to the chains. Along 
the chains this band is resolved into two sub-bands at 
3.08 and 3.26 eV. In addition, there is an extra band at 
2.65 eV, which is the lowest-energy excitation along the 
chains. 

The temperature evolution of the bands at 3.7 eV per- 
pendicular to the chains and at 2.65 and 3.45 eV along 
the chains can be also obtained from the dispersion anal- 
ysis. The total oscillator strength of the 3.7 (3.45) eV 
a bands is 10.76 (9.68) eV 2 . While the spectral weights 
of these bands are conserved and independent of temper- 
ature, their widths decrease gradually with decreasing 
temperature, as shown in the inset of Fig. 6(a). The 
dispersion analysis demonstrates a prominent strength- 



Table I: Parameters of Lorentz oscillators resulting from a dis- 
persion analysis in polarizations perpendicular and (along) 
the chains in NaCu202 measured at T = 7 K. £oo=2.00, 
(£00=1.77). 
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Sj(eV 2 ) 
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(2.65) 


(4.02) 


(0.60) 




3.37 (3.08) 
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0.57 (0.52) 
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Figure 7: Temperature-difference spectra Ae2{u) = 
£2(oj,T) — S2(w, 150 K) along the chains in (a) NaCu202 and 
(b) LiCuVC>4 at low energies. The red and blue arrows mark 
the positions of the two excitonic modes. 

cning of the optical band at 2.65 eV along the chains 
upon cooling below 150 K (see the inset of Fig. 6(b)), 
which is in full agreement with the measured dynamical 
scans shown in Fig 5(c). 

Close inspection of the temperature-difference spectra 
Aef (w) and Ae|(w) in the spectral region below 2.5 eV 
highlighted in Fig. 6(c) reveals a deviation of the fit in- 
corporating the bands listed in Table I (labeled as 'fit 1') 
from the experimental data indicated by stars. The devi- 
ation in both Aef (w) and Ae|(u;) spectra can be removed 
by introducing an additional band at 2.1 eV into the fit 
(labeled as 'fit 2') in Fig. 6(c). However, the weakness of 
the additional absorption band at 2.1 eV did not allow 
us to study its temperature dependence. 



III. DISCUSSION 

A. Optical response of the Cu 2+ 02 chains 

The optical response for both polarizations is com- 
posed of similar sets of oscillators, as revealed by the 
dispersion analysis. The only qualitative difference is the 
temperature-dependent band at 2.65 eV, which is clearly 
resolved at low temperatures along the chains. Compar- 
ison of its temperature evolution (Fig. 5(c) and inset of 
Fig. 5(b)) with magnetic susceptibility data^ points to 
the spin-controlled behavior of this band. When spin- 
correlations appear within the Cu 2+ C>2 chains below 150 
K, the magnetic susceptibility data deviate from Curie- 
Weiss behaviorJ^ Concurrently with the appearance of 
spin correlations, the intensity of the 2.65 eV band is 
enhanced, while it is virtually temperature independent 
above 150 K. 

The dispersion analysis at T = 7 K picked out a low- 
energy side band at 2.1 eV, see Fig. 6(c), apparent at 
7 K. Therefore, the lowest-energy excitation along the 
Cu 2+ C>2 chains has a doublet structure, as shown in Fig. 
7(a). Its characteristics including the energy positions, 
the partial spectral weights of the individual bands, and 



their evolution with temperature, are similar to those 
of the two-peak structure centered at 2.15 and 2.95 eV 
observed along the chains in LiCuVO*!^ This analogy 
indicates that these excitations have the same origin in 
both materials. Following arguments presented in Ref. 

we therefore ascribe the double-peak structure ob- 
served in NaCu202 to an exciton doublet formed by d 9 Lh 
and d 10 states generated by the NN (t x ) and NNN (t 2 ) 
hopping of electrons along the chains that emerges at 
U—V/l as a consequence of the long-range Coulomb inter- 
action between electrons on neighboring sites 1 = 1, 2^2 
The parameters characterizing the local Hubbard inter- 
action U = 3.2 eV and the long range Coulomb repulsion 
V = 1.1 eV can be obtained from the energies of these 
bands. Within this model the spectral weight for the 
first and second excitons are directly related to the cor- 
responding spin correlations^ 2 - 

N eff = -J^dlMSi ■ S i+l - 1/4), 1 = 1,2, (5) 

where di is the hopping length (d 2 = 2d\) and J; ~ 
4tf /(U — Vi). The large value of the antiferromagnetic 

(2) 

coupling constant, j r 2, leads to larger variation in 

as a function of temperature compared to N^X. The 
oscillator strengths of the individual bands in the differ- 
ence spectra of Ae|(w, 7 K) - Aejj(o>, 150 K) in Fig. 7 

implies a ratio N^l/N^fj « 24 in NaCu 2 2 and 11.5 
in LiCuVO.4. On the other hand, a ratio between the 
NN and NNN exchange integrals extracted form neu- 
tron scattering data is a ~ 5.5 in NaC^Oz-^, compared 
with the reduced ratio a ~ 2.7 in LiCuV04.— A precise 
quantitative determination of the exchange couplings in 
LiCuV04 is disputed3~— Nevertheless, the qualitative 
agreement between N^/N^ and a ratios for two com- 
pounds validates Eq. (5) for the spectral weight of the 
exciton bands. 

Specifically, the spin correlation functions in Eq.(5) are 
determined by the Heisenberg Hamiltonian of the under- 
lying frustrated J\-J 2 spin chainJ^ 

H s = Y; J ' £(3 ' Si+l - V4) + J[ 3 ' S i+ i, (6) 

l i i 

with the total NN exchange integral determined by the 
balance of two opposing contributions, J\ ot = J\ + j[ . 
Apart from the antiferromagnetic superexchange inte- 
grals J\ there is a substantial ferromagnetic coupling 
J-f that originates from a two hole excitation on the 
oxygen site that is not included in Eq.(5). This inter- 
action, however, plays an important role for the frus- 
trated magnetism of these spin chains. We may use 
now the additional experimental information of the rel- 
ative weights N^X/N^ff sa 24 at low temperature to- 
gether with the total exchange integrals obtained for 
NaCuaO^ J[ ot = -16 K and J|°* = J 2 = 90 K re- 
determine Ji and Jf separately. With help of an esti- 
mate of the correlation functions (Si-Si+i) ~ (—0.4) for 
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I = 1 (2) at zero temperature we find J\ ~ 39K, which 
implies for the ferromagnetic contribution j[ ~ — 55K. 



B. In-plane anisotropy and anomalous absorption 

The in-plane dielectric function of NaCu2C>2 measured 
with photon polarization parallel to the Cu 1+ planes ex- 
hibits strong absorption bands at 3.45 eV and at 3.7 eV 
perpendicular and parallel to the chains, respectively, 
which are unusual for strongly correlated electron sys- 
tems. The energy difference of 0.25 eV between their 
peak positions is a basic tendency also for other bands 
marked out by the dispersion analysis. The unique ex- 
ception is the 2.65 eV band, which is apparent only along 
the chains, as discussed in detail above. Since compounds 
with mixed-valent Cu atoms are scarce, these experimen- 
tal data are of special interest for studying the low-energy 
electronic excitations in cuprates. 

One of these mixed-valent compounds is the par- 
ent compound of a family of high-T c superconductors 
YBa 2 Cu 3 06, which shares some of its structural units 
with NaCu2C>2: Cu 2+ ions are centered in the CuC>4 pla- 
quettes forming the Cu02 conducting planes, and Cu 1+ 
ions construct the 0-Cu 1+ -0 dumbbells. Ellipsometric 
measurements with light polarized within the ab plane 
along the a-axis showed that the dielectric function of 
YBa2Cu3C>6 exhibits a sharp and intense peak at 4.1 
eV^ which disappears with increasing oxygen content^ 
Local density approximation (LDA) calculations assign 
the dominant absorption peak to intra-ionic transitions 
within Cu 1+ ions of the 0-Cu 1+ -0 complexes^ The ini- 
tial states were assigned to Cu 1+ 3d 3z 2_ 1 orbitals, and 
the final states consist of bonding combinations of Cu 1+ 
4p x orbitals and Ba 5d and 4/ orbitals. Although the 
Cu 1+ 3d 32 2_ 1 — > Ba transitions contribute only slightly 
to the optical matrix elements. Nonetheless, the analogy 
to the 0-Cu 1+ -0 transitions in NaCu2C>2 is not complete, 
because in this compound the oxygen atoms are shared 
with the Cu 2+ -02 chains. 

To our knowledge, the highest value of E2(u) among all 
transition- metal oxides has been observed in LiCu2 02, 
a compound that is isostructural and isoelectronic to 
NaCu2 02- In a recent ellipsometry study, Pisarev et al. 
reported that the real part of the dielectric permittivity 
of LiCu2C>2 exhibits an extremely strong and narrow ab- 
sorption peak at 3.27 eV. Overall, the in-plane dielectric 
response resembles our measured data and is composed of 
a similar set of Lorentz oscillators. However, Li-Cu chem- 
ical substitution, a 10% nonmagnetic LiCuO impurity 
phase and the twinned nature of LiCu2C>2 crystal o 26 i 27 
obliterate the dielectric anisotropy in the ab plane. Re- 
viewing the optical data of the large family of cuprates 
with 0-Cu 1+ -0 complexes, Pisarev et al~- pointed out 
a relationship between the Cu 1+ -Cu 1+ and Cu 1+ -0 2 ~ 
bond lengths and the position and intensity of the sharp 
peak, which implies that this feature originates from the 
dumbbells. Arguing on the basis of the peculiar inten- 
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Figure 8: (a) Energy band dispersion calculated by the spin- 
restricted LSDA approximation for NaCu2C>2 with partially 
filled Cu 2+ 3d (red) and unoccupied Cu 1+ hybrid 4p and 
3d 3z 2_! (blue) states, (b) DOS projected on the Cu 2+ 3d 
(red) and Cu + 3d (blue) atomic states of NaCu2 02. The 
Fermi level is at zero energy. 
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Figure 9: Calculated diagonal elements of the (a) erf (uj) and 
(b) <j^(uj) optical conductivity of NaCu202 (black solid line) 
and its decomposition into transitions to the Cu 2+ (red dash- 
dotted line) and the Cu 1+ (blue dashed line) final states. 



sity of the 3.27 eV band, the authors proposed an exciton 
model and attributed the observed anomaly to a strong 
Cu 1+ 4p-3<i electron-hole interaction along with a strong 
crystal-held splitting of the excited states. 

In order to interpret the observed anisotropy and clar- 
ify the origin of the intensive absorption peak, we per- 
formed band structure calculations of NaCu 2 02 along the 
high symmetry directions of the Brillouin zone. 
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Photon energy (eV) 

Figure 10: (a) Partial contribution of transitions from the occupied Cu 1+ 3d 3z 2_ 1 to the unoccupied Cu 1+ 4p x (blue dash- 
dotted line) and Ap y (red dashed line) states to erf (uj) and a^(ui) within spectral range 2.35—3.5 eV. (b) The belts of transitions 
from Cu 1+ 3d z 2_ 1 to 4p x (blue) and to 4p y (red) states, partially contributing af(uj) and <Ji(uj) in (a) and the energy band 
dispersion in fat-band representation. The size of black squares, blue and red circles in (b) is proportional to the partial weights 
of the Cu 1+ 3d z 2_ 1 , 4p x and 4p y states in the Bloch wave function, respectively. 



C. Band Structure Calculations 

Band structure calculations were performed using a 
linear-muffin-tin orbital method in the atomic sphere 
approximation 2 ^ within the local spin density approxima- 
tion (LSDA) starting from the known crystal structure 
of NaCu 2 2 . 

The results of the spin-restricted LSDA calculations 
along with the resulting partial electronic densities of 
states (DOSs) are shown in Fig. 8. The Cu 1+ and Cu 2+ 
3d DOS are shown in Fig. 8(b), as they are the main 
contributors to the total DOS. Since the LSDA results in 
a metallic state, narrow partially filled bands correspond- 
ing to Cu 2+ 3d orbitals with xy, xz symmetry cross the 
Fermi level and give a DOS peak at Ep. The valence 
band between —4 eV and —1 eV below the Fermi energy 
Ep is a mixture of Cu 1+ and Cu 2+ 3d states of different 
symmetries strongly hybridized with O 2p states. The 
modestly intense band at near 2 eV above the Fermi en- 
ergy shown by blue lines in Fig. 8(a) is mainly formed 
by unoccupied Ap and 3d 3z 2_ 1 states of Cu 1+ . 

The calculated optical conductivity spectra a\ (w) and 
(Ji(uj) obtained from the band structure are shown in 
Fig. 9. The sums of all the interband transition con- 
tributions to the optical conductivity (the total optical 
conductivity) are shown by black solid lines. Blue dashed 
and red dash-dotted lines indicate the spectral weights 
associated with transitions to the unoccupied Cu 1+ and 



Cu 2+ states, respectively. The calculated optical re- 
sponse below 2.4 eV originates entirely from transitions 
involving Cu 2+ final states, as discussed in detail below. 

The calculated spectra above 2.4 eV, contributed 
mostly by transitions into the Cu 1+ states, conform to 
the experimentally observed response in the 7 region, 
which exhibits two dominant bands at 3.45 and 3.7 eV 
(a region Fig. 3) on top a broad background (/? and 7 
regions). While the dispersion analysis of the experimen- 
tal data separates the 7 zone into bands with half-widths 
~ 0.5 eV (Fig. 6(a) and 6(b)), the LSDA gives sharper 
structures. However, within this spectral range, the aver- 
age calculated optical conductivity ~ 2.5xl0 3 17 _1 cm~ 1 
is in agreement with the measured data. The intensities 
of both calculated and experimental spectra are reduced 
above 5 eV. 

To find out the origin of the observed in-plane 
anisotropy, appearing as a shift of the spectra by 0.25 
eV, we analyzed the interband transitions contributing 
to the optical conductivity above 2.4 eV. The decom- 
position of the calculated cr^(u>) and <j\{uj) spectra into 
separate transitions revealed that the Cu 1+ 3d 3z 2_ 1 or- 
bitals serve as initial states for transitions within the 
spectral range 2.4—3.5 eV in Fig. 10(a). The final states 
involved are of Cu 1+ 4p character with x and y symme- 
tries. The respective transitions and dispersions of bands 
in a fat-band representation are illustrated in Fig. 10(b). 

As shown in Fig. 10(a), the erf (a;) spectrum around 3 
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Photon energy (eV) 

Figure 11: (a) Partial contribution of transitions from the occupied Cu + 3d xz (blue dash-dotted line) and and 3d yz (red 
dashed line) to the unoccupied Cu 2+ 3d xy states to <t|(w) and within the spectral range from 1—2 eV. (b) The belts of 

transitions from Cu 1+ 3d xz (blue) and Cu 1+ 3d yz (red) to Cu 2+ 3d xy states, partially contributing to o\{lS) and <Ji(lu) in (a) 
and the energy band dispersion in fat-band representation. The size of black squares, blue and red circles is proportional to 
the partial weights of the Cu 2+ 3d xy , Cu 1+ 3d xz and Cu 1+ 3d yz states in the Bloch wave function, respectively. 



cV is, mainly, formed by transitions into the Cu 1+ 4p x 
states (blue dash-dotted lines) , while transitions into the 
4p y states (red dashed lines) dominantly form the <t^(uj) 
spectra. Because of the large dispersion of these states, 
the transitions give a rather flat shape of the optical con- 
ductivity in this spectral range and spread out up to 5 
cV. Above 3.2 cV the calculated spectra arc formed by 
the Cu 1+ 3d unoccupied states and O 2p — > Cu 2+ 3d in- 
terband transitions, contributing equally to both a^(u>) 
and erf (a;). 

Let us turn now to the interband transitions associ- 
ated with the Cu 2+ 3c? final states, giving, in particular, 
the sharp peaks near 1.5 eV in Fig. 9. The decompo- 
sition into partial contributions reveals that these peaks 
are formed by transitions from the occupied Cu 1+ 3d XZtVZ 
bands to the Cu 2+ 3d states with xy symmetry, see Fig. 
11(a). The high intensity and sharpness of the peaks are 
explained by the low dispersion of the Cu 2+ 3d xy states, 
lying at the Fermi level in the LSDA calculations, and the 
Cu 1+ 3d xz:yz bands near the T points of the Brillouin- 
zone in Fig. 11(b). The initial states with different sym- 
metries, xz and yz, contribute differently to the optical 
conductivities trj L (w) and c|(u;), as indicated in Fig. 11. 
While these transitions have predominantly d character, 
the contribution of p states to the Cu 1+ 3d XZtVZ bands 
is sufficient to provide a significant transition probability 
through optical dipole matrix elements. 

The above analysis of the interband transitions, which 



are responsible for the strong anisotropy of the optical 
properties of NaCu2C>2, is based on the spin-restricted 
LSDA calculation which gives a metallic solution with 
four Cu 2+ dry-derived bands crossing the Fermi level 
(Fig. 8). Although a minute gap of less than 0.1 eV 
opens in LSDA calculations for spin-spirals (not shown), 
it is still much smaller than the experimental gap of 2 
eV. The reason for the discrepancy is that the strength 
of electronic correlations in the Cu 2+ 3d shell is strongly 
underestimated within the LSDA. When properly ac- 
counted for, the on-site Coulomb repulsion Udd would 
split the half-filled Cu 2+ d xy bands into occupied lower 
and unoccupied upper Hubbard bands and open an in- 
sulating gap. Then, the sharp peaks of the calculated 
optical conductivity, which appear at ~1.5 eV due to the 
interband transitions involving the Cu 2+ d xy final states 
(Fig. 8), would shift to higher photon energies improving 
the agreement with the measured spectra. 

In order to illustrate the effect of the electronic corre- 
lations on the optical conductivity of NaCu202 we recal- 
culated its band structure and optical spectra using the 
LSDA+C/ method^. The LSDA+C/ calculations were 
performed assuming ferromagnetic order of the Cu 2+ mo- 
ments. The value of the exchange integral J was fixed to 
the LSDA value of 1 eV, and the on-site Coulomb repul- 
sion Udd was varied from 2 to 8 eV. Before presenting the 
results of the LSDA+C/ calculations it is worth recalling 
that when the so-called atomic limit is used for the dou- 
ble counting terrn^ and non-spherical contributions to 
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Udd and J are neglected, the expression for the orbital 
dependent LSDA+J7 potential V ail which is to be added 
to the LSDA potential, becomes particularly simple: 

V ai = {Udd - J) Q - , (7) 

where n a i is the occupation of i-th localized orbital with 
the spin a. One immediately notices that the main effect 
of LSDA+t/ is to split occupied {n a i ~ 1) and unoccu- 
pied {n a i w 0) states by shifting the former by {Udd—J)/2 
downwards and the latter by the same amount upwards 
with respect to their LSDA energy position. 

Since all the Cu 2+ d states, except for d xy ones, 
are completely filled already in LSDA, they shift by 
{Udd — J)/2 to lower energies when Udd is applied. The 
half-filled d xy states contribute substantially to bonding 
Op-Cu d states near the bottom of the valence band be- 
cause of strong (T-type hybridization with O {p x ±p y )-\/2 
states. In LSDA+J7 calculations the majority spin d xy 
state becomes fully occupied and moves by {Udd — *0/2 
downwards. The occupation of the minority spin d xy 
state, however, does not go to zero, as one would expect 
for a formally unoccupied state, but instead remains close 
to 0.5 due to the strong d xy contribution to the occupied 
O p— Cu d bonding states. According to Eq. (7), the 
LSDA+J7 potential acting on the minority spin Cu 2+ d xy 
state is much less than {Udd — J)/2, and the energy of the 
bands formed by these states remains close to the LSDA 
value. In order to avoid such asymmetric splitting of 
the d xy bands, we performed the LSDA+[7 calculations 
with the occupation numbers of the majority and minor- 




Photon energy (eV) 

Figure 12: Optical conductivities crj'(w) and (b) Oi{uj) cal- 
culated with Udd=6 eV acting (a) on the Cu 2+ and Cu 1+ d 
states and (b) on the Cu 2+ d states only. 



ity spin Cu 2+ d xy states fixed to 1 and 0, respectively. In 
this way the corresponding majority and minority spin 
bands arc shifted by ±{Udd — </)/2 with respect to their 
LSDA position. 

The optical conductivities cr^{io) and o\(uS) calculated 
with U^ d u =6 eV acting on the Cu 2+ d states are pre- 
sented in Fig. 12 (b). With this value of U%/ + , the 
unoccupied bands formed by the minority spin Cu 2+ d xy 
states lie ~2 eV above the Fermi level and, in contrast 
to the LSDA result, the LSDA+J7 band structure is in- 
sulating with a gap of about 1.8 eV. The sharp absorp- 
tion peaks, which are caused by inter-band transitions to 
the final states of the Cu 2+ d xy character and appear at 
~1.5 eV in the LSDA spectra (Fig. 9), shift to 3 eV in 
much better agreement with the experiment. The peaks 
show strong polarization dependence, with the <j\ peak 
being 0.3 eV lower than the erf one. When Ug u of 
6 eV is applied also to the d states of Cu 1+ ions, the 
bands formed by Cu 1+ d yz and d zx states, which are 
the initial states for the interband transitions responsi- 
ble for the sharp peaks, shift to lower energies. In this 
case all n a i ~ 1, and LSDA+J7 mimics to some extent 
the effect of the so-called self-interaction corrections^! 
by shifting all the Cu 1+ d states downwards. This brings 
the conductivity peaks even closer to their experimen- 
tal positions. Nevertheless, the magnitude of the cal- 
culated optical conductivity remains lower than in the 
experiment. One would also expect some narrowing of 
the unoccupied Cu 2+ d xy bands if the LSDA+C/ calcu- 
lation were performed for the incommensurate magnetic 
structure observed experimentally! 14 i 1 5 This could result 
in sharpening of the peaks and an increase of the calcu- 
lated conductivity. 

Our calculations suggest that the origin of the anoma- 
lously strong absorption peaks observed at 3.45 and 3.7 
eV can be ascribed to the interband Cu 1+ 3d —> Cu 2+ 
3d transitions. The observed shift in the peak positions 
between two polarizations arises from different contribu- 
tions of transitions from the initial d zx and d yz states 
to the optical conductivity, as shown in Fig. 11. The 
experimentally observed anisotropy of the background, 
assigned to the Cu 1+ 3d — > Ap transitions, is due to the 
different contributions of transitions to the final states 
with p x and p y characters to the optical conductivities 

<Ji{lo) and cr|(a;), as shown in Fig. 10. 

IV. CONCLUSIONS 

In our spectroscopic ellipsometry study of NaCu202, 
we observed that the spectra of the dielectric function 
for light polarized parallel to the Cu 1+ planes exhibit a 
strong in-plane anisotropy of the interband excitations. 
The absorption edge for polarization along the Cu 2+ 02 
chains is formed by a two-peak structure centered at 
2.15 and 2.65 eV. This feature is absent in the other 
polarization and bears striking resemblance to the one 
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observed in the single- valent Cu 2+ C>2 chain compound 
LiCuV04^ Our findings suggest that these modes are 
a generic feature of Mott-Hubbard insulators with edge- 
sharing copper-oxide chains. Based on theoretical con- 
siderations, we have also shown that an exciton doublet 
is expected to emerge as a consequence of the long-range 
Coulomb interaction at energies U—V and U — V/2 upon 
cooling into a temperature range where substantial spin 
correlations become established. Identification of this 
exciton doublet with the experimentally observed two- 
peak feature allowed us to determine the local Hubbard 
interaction parameter U = 3.2 eV and the long range 
Coulomb repulsion parameter V = 1.1 eV. 

The quantitative information about elementary excita- 
tions in insulating chain cuprates gained from the present 
study may deepen our understanding of the electronic 
structure and phase behavior of doped cuprates as well. 
In the doped Cu02-planes of the cuprate superconduc- 
tors the band width is significantly larger, and hence exci- 
tons are not expected to appear. Nonetheless, the quanti- 
tative description of the effect of the long-range Coulomb 
interaction obtained here may facilitate realistic calcula- 
tions of the properties of charge density wave and striped 
states in doped compounds with one-dimensional^ and 
two-dimensional electronic structure. 



We have also observed strong and sharp absorption 
bands peaked at 3.45 eV and 3.7 eV for light polarization 
along and perpendicular to the Cu 2+ C>2 chains, respec- 
tively, dominating the spectra and superimposed on a fiat 
and featureless plateau background above the absorption 
edge. Based on density functional calculations, we con- 
clude that the major contribution to the background re- 
sponse comes from the intra-atomic Cu 1+ 3c? — > 4p tran- 
sitions within the 0-Cu 1+ -0 dumbbells. The experimen- 
tally observed anisotropy of the background is explained 
due to the different characters of the final Ap x and Ap y 
states. Pisarev et al. have attributed this anomalous 
absorption feature to exciton formation of the Cu 1+ Ap 
electron - 3d hole pairsJ^ We propose an alternative ex- 
planation based on consideration of electron correlations 
in the LSDA calculations. Our results indicate that the 
anomalous peak superimposed on the background ab- 
sorption can be assigned to transitions between bands 
formed by Cv} + 3d xz /d yz and Cu 2+ 3d xy orbitals, which 
arc strongly hybridized with O-p states. In this approach, 
the observed ~ 0.25 eV shift in the peak positions be- 
tween two polarizations arises from anisotropic contribu- 
tions of transitions from the initial Cu 1+ 3d xz and 3d yz 
states to the optical conductivity. 
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